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Physics of the interior of a black hole with an exotic scalar matter 
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ABSTRACT We use a numerical code to consider the nonlinear processes arising when a 
Reissner-Nordstrom black hole is irradiated by an exotic scalar field (modelled as a free mass- 
less scalar field with an opposite sign for its energy-momentum tensor). These processes are quite 
different from the processes arising in the case of the same black hole being irradiated by a pulse of 
a normal scalar field. In our case, we did not observe the creation of a spacelike strong singularity 
in the T-region of the space-time. We investigate the antifocusing effects in the gravity field of the 
exotic scalar field with the negative energy density and the evolution of the mass function. We 
demonstrate the process of vanishing of the black hole when it is irradiated by a strong pulse of an 
exotic scalar field. 

PACS numbers: 04.70.Bw, 04.25.D-, 95.36.-|-x 



I. INTRODUCTION 

The internal structure and physics of black holes (BH) 
has been the subject of researches during many years 
W\ . A powerful tool for these investigations is to consider 
a spherical, charged, nonrotating BH which is nonlinearly 
perturbed by a selfgravitating scalar field. This toy BH 
model is not very realistic but it shares many proper- 
ties, including causal structure, with the more realistic 
rotating BHs. 

This toy model has been used in the paper to an- 
alyze the physics of the interior of a BH in the case of 
irradiation by a normal massless scalar field. 

The observational apparent acceleration of the uni- 
verse suggests the presence of a matter field which vi- 
olates at least the strong and perhaps also the weak en- 
ergy condition and one of the simplest examples of such 
energy-conditions-violating matter is a free scalar field 
with a reversed sign of its energy-momentum tensor, we 
call this an exotic scalar field. Such an exotic scalar field 
has previously been used in many works for construc- 
tion of wormholes in cosmological models with nontrivial 
topologies, see for example [SJ,!!!,!!^- It would thus be 
interesting to explore possible implications of such a field 
to black hole physics. 

The goal of this paper is to perform such an analysis, 
where we will model the exotic matter as described above. 
We will sec that the physics of the interior of a BH which 
is nonlinearly perturbed by such an exotic scalar field is 
quite different from the physics of a BH perturbed by 
normal scalar field. 



Rather recently, some aspects of the problem of the 
internal structure of a BH with unusual scalar fields have 
been discussed in the papers [40l - |4^ . see also references 
therein. 

In the paper ^i^l the authors discussed the interior 
structure of a BH arising as a result of the collapse of a 
spherical charged scalar field shell, taking into account 
the Hawking radiation. The problem was treated numer- 
ically. 

In the paper [4l[ some problems of the inner structure 
of a BH were discussed using an analytical approach. 

Finally, in the paper p^ . the authors explored nu- 
merically the accretion of a phantom scalar field with a 
nonzero potential onto a noncharged BH. Especially the 
variation of the BH area were investigated. 

Our approach is different from the above papers. We 
will address the following points: 

• What physical processes arise in the case of irradia- 
tion of a charged BH by an exotic scalar field and what 
is the difference from the case of irradiation by a normal 
scalar field. 

• What kind of singularities arise inside of a charged 
BH in the case of irradiation by an exotic scalar field. 

• How is the the charged BH transformed into an ob- 
ject without horizons in the case of the BH being irradi- 
ated by stronger pulses of exotic scalar radiation? 

• What are the proper Penrose diagrams for both sur- 
viving and nonsurviving (charged) BHs in the case of 
irradiation by an exotic scalar field. 

We will study the nonlinear processes inside the BH us- 
ing a stable, second order accurate numerical code with 
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adaptive mesh refinement capabilities. This code was 
described and tested in details in [S^l, see also the Ap- 
pendix El 

This paper is organized as follows: 

In section [III Einstein's equations are written for a 
spherical BH with a fixed magnetic charge and massless 
scalar fields, both normal and exotic. In section Hill we 
formulate the initial values for our computations. In sec- 
tion |lVl we discuss the mass function and some nonlinear 
processes inside a BH. In section |Vl we discuss the prob- 
lem of the origin of strong singularities in a BH. In sec- 
tion lVIl we describe our numerical model. In section lVlIi 
we remind of the physical processes inside a charged BH 
irradiated by a pulse of the normal scalar field. In sec- 
tion I Villi we analyze the results of our computations 
which describe the physics of the interior of a BH with a 
magnetic field irradiated by a pulse of exotic scalar radi- 
ation. In section HXl we discuss the case of irradiation of 
the BH by both normal and exotic radiation pulses and 
we summarize our conclusions in section |Xl Finally, in 
appendix A, we presents results of convergence tests of 
the numerical code. 



II. FIELD EQUATIONS FOR THE SPHERICAL 
MODEL 



We wish to investigate the evolution of a spherical BH 
with a fixed magnetic or electric charge q (i.e. Reissner- 
Nordstrom metric) under the action of pulses of selfgravi- 
tating, massless scalar fields, both normal $ (with a posi- 
tive energy density, e > 0) and exotic ^ (with a negative 
energy density, e < 0). The general equations for the 
analysis were written in (s^ ]. 

The line element in double-null coordinates can be 
written as 

ds^ = ~2e^''^''''''>dudv + r^{u,v)dn^, (1) 

where cr(u, v) and r(u, v) are functions of the null coor- 
dinates u and V (in- and out-going respectively). The 
nonzero components of the Einstein-tensor are: 



electro-magnetic field, i.e. : T^^ = T*^ + T*^ + T^™ 
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The energy-momentum tensor can be written as a sum 
of contributions from the exotic scalar field ^ (with the 
negative energy density), from the normal scalar field $ 
(with the positive energy density) and from the ordinary 
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The u — u. V — V, u — V and 9 — 9 components of the 
Einstein equations (with c — 1, G — 1) respectively are: 
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The scalar fields satisfy the Gordon-Klein equation 
V^V^* = and V^V^* = 0, which in the metric (P) 
become: 

+ - (?-,^*.u + r„^'.„) = (10) 
r 

^,uv + - {r,v^.u + r „$,„) = (11) 
r 

Equations (|8l)- (|lip are evolution equations which are sup- 
plemented by the two constraint equations © and (O. 
It is noted that none of these equations depend directly 
on the scalar fields ^ and $ but only on their derivatives, 
i.e. the derivative of the scalar field is a physical quantity, 
while the absolute value of the scalar field itself is not. 
Specifically we note the T*„ = -(5',„)V(47r) and T* = 
— (5* /{All) components of the energy-momentum ten- 
sor @ and T*„ = ($,„)V(4^) and T* = (*,.)V(4^) 
components of the energy-momentum tensor ^ , which 
are part of the constraint equations. Physically r„„ and 
Tyy represents the fiux of the scalar field through a sur- 
face of constant v and u respectively. 



III. INITIAL- VALUE PROBLEM 

We wish to evolve the unknown functions r{u,v), 
a{u,v), $(u,ti) and '^(u,v) throughout some computa- 
tional domain. We do this by following the approach 
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of [11 [13, m H, m H, m 113, S, m to numeri- 
cally integrate the four evolution equations ((8|)- (fTT|) along 
the characteristics. These equations form a well-posed 
initial-value problem and in the double-null coordinate 
we use, it is most natural to use the characteristic initial- 
value formulation in which the initial values of the un- 
known functions are specified on two initial null seg- 
ments, namely an ingoing (v ^ vq = constant) and an 
outgoing {u = uq — constant) segment |43j . 

Our code is a free evolution code, i.e. we impose 
the constraint equations © and ([7]) only on the ini- 
tial segments, consistency of the evolving fields with the 
constraint equations is then ensured via the contracted 
Bianchi identities [43|. However we use the constraint 
equations throughout the domain of integration to check 
the accuracy of the numerical simulation. 

We are free to choose the distribution of the $ and 
^'-fields on the initial null segments, our choices are de- 
scribed in section lVTl We are also free to choose r{u, v) on 
the initial surfaces, this merely expresses the gauge free- 
dom associated with the transformation u — > u{u),v — > 
v{v) (the line element (dJ and the equations (p])- pT|) are 
invariant to such a transformation). Hence, the only 
variable left for us to determine on the initial surfaces 
is a. This can easily be found by integrating the con- 
straint equations eq. ^ and ([7]), which ensures that the 
constraint equations are satisfied on the initial hypersur- 
faces. Specifically on the outgoing u — uq hypersurface 
it is found by the integral: 

( \ i \^ } r...-r{^^,)^ + r ($,,)^ ^ 
a[uo,v) ^ a[uo,vo) + / dv 

J ^ 

vO 

where cr(uo, vq) is an integration constant. We follow [27j, 

m, m, m El and set cr(uo,wo) = 1^(71) throughout 
the paper. The initial data along the ingoing surface, is 
set in a similar manner. 

Hence, by specifying a distribution of the scalar fields 
$ and ^I* on the initial null segments, choosing a gauge, 
integration constant a{uQ^vo), charge and mass parame- 
ters "g" and "m" (these relate to the gauge, see eq. (fT3|) 
and section lVll) we can specify complete initial conditions 
on the initial null segments. Using the numerical code 
from ['s^l, we can then use the evolution equations, eqs. 
(I5))- ([TT|1 . to evolve the unknown functions along the char- 
acteristics, throughout the computational domain. Our 
specific choices for initial conditions are described in sec. 

ED 

Finally, we wish to emphatize that a consequence of 
evolving along the characteristics in our double-null co- 
ordinate, spherical model, is that the specification of ini- 
tial data along in- and outgoing null surfaces fully defines 
the initial-value problem. Thus, there is no need to spec- 
ify any additional boundary conditions, the initial values 
along the null surfaces are the boundary conditions. For 
further details on the numerical code, see [sj] and refer- 
ences therein. 



IV. MASS FUNCTION 

The evolution of the interior of a BH with a scalar field 
is highly nonlinear. An important characteristic of this 
evolution is the mass function. This function represents 
the total mass (without the ma gne tic field) in a sphere 
of radius r{u,v) (see [13, S IMSl ) ■ In the metric ([J), 
the mass function is (see [23,1301): 

In the metric 

ds"^ = gtt dt^ + grr dr'^ + dil'^ (14) 
the mass function m has the form 

(1 + ^-5.-/) (15) 

There are two physical processes which can lead to a 
nonlinear change of the mass function: 

1. The mass m inside a sphere can change because of 
the work of the pressure forces acting on the surface of 
the sphere. 

2. Mass inflation can occur [Toj . 

Both processes have been described in [H. We want 
to emphasize that in the case of the exotic scalar field 
we should take into account the negative energy density 
and negative radial pressure of the exotic scalar field. 
Of course, in addition to the above mentioned processes, 
there is also the trivial process of matter fiowing into 
the sphere. We will discuss the properties of the mass 
function in our model in the subsequent sections. 

One more important nonlinear process is the focusing 
effect caused by the gravity of pulses of opposite fiuxes 
of radiation (see [ll). In our case of radiation pulses 
with negative energy densities, it is rather a defocussing 
effect. Its manifestation will be discussed below. 



V. STRONG SINGULARITY 

In the case of nonlinear perturbations of the initial BH 
by the normal $ scalar field, a strong spacelike singularity 
arises in a T-region as a result of the nonlinear processes 
(see for example [H, [H, i4^] ) . 

We want to remind of the classification of R- and T- 
regions in space-times with spherical symmetry (46l - l48| . 
The i?-regions are regions where world-lines r = const are 
timelike. The T+-regions are regions where world-lines 
of r = const are spacelike and the direction to larger 
r is in the direction to the future. The T_-regions are 
regions where world-lines of r = const are spacelike and 
the direction to smaller r is in the direction to the future. 
The R- and T-regions are separated by apparent horizons 
which are defined by locations where 4^ = or = 0. 

^ av au 

The strong singularity under discussion exists together 
with a weak null singularity (instead of the inner horizon 
of a Reissner-Nordstrom BH 
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In the case of nonlinear perturbations of the initial BH 
by the exotic vl/-ficld the situation is quite different. 

First of all, the analysis in [4^ have demonstrated that 
in the case of a VE* field in the vicinity of the singularity 
r = 0, we can neglect the influence of the exotic scalar 
field 5* comparative to the infiuence of the magnetic field. 
This means that the singularity r = is of a Reissner- 
Nordstrom type. It is a timelike singularity and it is in 
an R-region. A spacelike strong singularity r = in a 
T_-region is not possible. 

Secondly, we will see in the subsequent sections that 
under the action of the irradiation by the exotic scalar 
field, the T_-region becomes smaller and smaller when 
the intensity of the irradiation increases, even in the case 
when the BH still survives. We can see practically all 
of the T^-region in our computational domain and def- 
initely there are not any strong singularities in the r_- 
region. On the other hand, a strong timelike singular- 
ity of the Reissner-Nordstrom type must exist in the R- 
region outside of our computational domain. 

We will consider the problems of the singularities in 
some detail in the following sections. 

Finally, we wish to note that no special treatment is 
needed in order for our numerical code to successfully 
handle the spacelike r = singularities that occur in our 
computations. The characteristic evolution of our nu- 
merical code means that only points which are causally 
connected to a singularity can "see" the singularity, but 
since no points are causally connected to a spacelike r = 
singularity, they pose no problem. Timelike and coordi- 
nate (r — 0) singularities, on the other hand, do pose 
a problem in principle as these points are causally con- 
nected to the rest of the space-time. In practice, however, 
it is not a problem since we choose our initial conditions 
in such a way that they do not include the coordinate 
singularity (cf. sec. IVip and any timelike singularity are 
outside of our computational domain. 

For further details on the workings of our numerical 
code, please see [13] and references therein. 



VI. THE MODEL 

In the subsequent sections we investigate the fully non- 
linear processes in the spherical charged BH irradiated by 
a pulse of the exotic scalar field ^ (and in section IIXI ir- 
radiated by both 5' and $ pulses). For that we use the 
results of our numerical simulations. Our numerical code 
is thoroughly described and tested in [sJl, here we just 
mention that the code uses Adaptive Mesh Refinement 
(AMR) in both u— and w— coordinates and is second or- 
der accurate. About the convergence of the code in the 
case of the present model, see Appendix [X] 

Our choice of the initial values corresponds to the fol- 
lowing physical situation; There is a BH with a magnetic 
field and at some distance from the horizon, at the ini- 
tial moment, there is a rather narrow spherical layer of 
an in-falling scalar field. 



Throughout this paper, prior to the irradiation by the 
scalar pulses, the BH has initial mass mo = 1 and charge 
q — 0.95mo. For the simulations, the domain of integra- 
tion as a rule is 5 < u < 30, < m < 30 and we choose 
a standard gauge in which r is linear in v and u on the 
initial null segments. Specifically, our gauge choice is: 

r{uQ,v) = v r{u,Vo) ^ r{uo,Vo) +r^uU (16) 

where r.„ is determined from the mass function, eq. (jl3l) . 
at the initial moment uo,vq, to ensure consistency. We 
note that this computational domain does not include 
the r — coordinate singularity. 

For our numerical experiments we specify various ini- 
tial conditions. Along the ingoing border vq — 5, the ini- 
tial values for all simulations correspond to the (Reissner- 
Nordstrom) BH solution with mo = 1 and q = 0.95mo 
without any additional fields. Along the outgoing border 
uq = 0, we vary the initial conditions in such a way as 
to imitate some physical fiuxes of the fields $ and ^E* into 
the charged BH. 

We specify the pulses of the scalar fields in the form: 

<^Ano, v) = sin^ (n ^ " ) , (17) 

"f.viuo, v) = sin^ [tt-^^—^^] , (18) 
\ - i;*o / 

where Aq, and A>p, v^o, v^q, w>pi are constants. 
Outside the interval (f$Oj v^i) we have ^,y(uo, v) = 0. 
Outside the interval (w^Oi^^ii) we have ^',„(uo,f) = 0. 

If the fluxes through the two initial surfaces uo and 
vq have been chosen, as was described in section ITTTl all 
other initial conditions are determined by our choice of 
gauge and the constraint equations. 

VII. PHYSICAL PROCESSES INSIDE OF THE 
CHARGED BLACK HOLE IRRADIATED BY A 
PULSE OF NORMAL RADIATION 

We start from reminding of what happens if we irradi- 
ate a charged BH with a pulse of the normal scalar field 
with a positive energy density [H, [H, 113, ^2^, 'ss'] . 

We described the corresponding processes in detail 
in [33], here we summarize main effects: 

1) When the positive energy pulse crosses the outer 
apparent horizon (OAH), the horizon becomes bigger. 
When the pulse crosses the inner apparent horizon 
(lAHl), it becomes smaller. 

After the passage of the pulse, the mass function in- 
creases by an amount corresponding to the energy con- 
tained within the pulse. These effects are the trivial con- 
sequences of mass being pumped into the BH by the ini- 
tial pulse. 

2) Next is the focusing effect. Any test photons prop- 
agating along the null-geodesic lines u = const and 
1; = const are under action of the gravity of the fluxes Tyy 
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and Tuu respectively. The last arises as a result of the 
scattering of the ingoing Tyy radiation by the space-time 
curvature. In the subsequent evolution these T^^ and T„„ 
fluxes are converted into one another due to the curva- 
ture of the space-time. The gravity of Tyy and Tuu leads 
to a focusing effect on the test photons. In the absence of 
the scalar radiation, outgoing photons along u = const 
slightly inside lAHl will go to greater r as w increases 
in the Reissner-Nordstrom solution. With the existence 
of the scalar radiation, a similar outgoing ray will now, 
because of the focusing effect, go to smaller r and gener- 
ate a maximum 4^ I ^ = in the u — v coordinates. 

av I r— const 

This corresponds to the position of the apparent horizon. 
As a result the position of lAHl changes dramatically. It 
goes up sharply with increasing v coordinate. 

3) Now we consider the properties of the mass function. 
For late (big) u the mass function along the lines u — 
const demonstrates a sharp increase with growing v. This 
is the result of the compression which tends to infinity. 
As a consequence, a strong singularity r = arises and 
these lines, u — const, come to the singularity r = 0. The 
Kretschmann scalar {K — Rap-^sR"'^^^), which is the 
curvature characteristic of the space-time, also increases 
sharply along these lines, due to the infinite compression. 

For smaller u the lines u = const go to the Cauchy 
horizon lAHl at ?; — > oo. The increase of the mass func- 
tion along these lines is mainly due to the mass inflation 
and is exponential when v oo. The mass inflation is 
the basic process in the formation of the weak singularity 
instead of the inner smooth Cauchy horizon. 

The weak null singularity shrinks down due to the fo- 
cusing effect from the gravity of the Tuu radiation and 
finally comes down to r = 0. 



VIII. PHYSICAL PROCESSES INSIDE OF A 
CHARGED BLACK HOLE IRRADIATED BY A 
PULSE OF EXOTIC SCALAR RADIATION 

Now let us compare the processes described in the 
previous section with the processes in the case when a 
Reissner-Nordstrom BH is irradiated by an exotic scalar 
field. 



A. The case of survival of the BH 

We start from the case when the power of the exotic 
radiation pulse is rather weak and the BH survives after 
being irradiated. Let us first consider the following case: 



0,; Axt^O.Ol (wvKo 



i;*i=9). (19) 



We will see that total amount of negative mass pumped 
into the BH in this case is rather small, it survives but 
the positions of the apparent horizons change. 

It is clearly seen from fig. 1(a) and fig. 1(b) that the 



result of the processes with the negative energy density 
pulse is opposite to the case of a positive energy density 



pulse described above. When the pulse of the negative 
energy density crosses the OAH, the horizon becomes 
smaller and on the corresponding plot (Fig. 1(a) I it can 



be seen to be going to higher u (corresponding to smaller 
r, cf. eq. (dH)) in the region 5 < w < 9. 

Before the pulse at w < 5, the mass of the BH was 
m = 1 and the position of the OAH corresponded to 
u = 23.190. After the end of the process, the final posi- 
tion of the OAH corresponds to m w 23.206 (fig. 
and the final asymptotic external mass is m w 0.998 
we will see in fig 



1(a)) 



as 



2(a)). 



We formulated the initial conditions along the outgo- 
ing null surface u = 0. When the ingoing radiation pulse 
propagates from w = to the OAH, it experiences scat- 
tering because of the curvature of the space-time, thus 
being converted into outgoing radiation. At later times, 
that scattered radiation is then rescattered and being 
converted back into ingoing radiation and so on and so 
forth. The tail of this scattered radiation thus reaches 
the OAH at later times. This can be observed in fig. 



1(a) I where the OAH can be seen to change its position 



in two steps. 

We would like to emphasize that the scattering and the 
rescattering are the result of interaction of the propagat- 
ing field with the curvature of the space-time background 
but not any boundary conditions. 

The first step, when the initial pulse passes into the 
BH, can be observed between 5 < < 9, and the second 
step (caused by the scattered radiation) can be observed 
between 9 < u < 11. 

When the pulse crosses the lAHl, the effects are again 
opposite to those of an energy pulse with a positive en- 
ergy density and the horizon thus becomes bigger (goes 
to smaller u) in the region 5 < w < 9 (Fig. 

In the region w > 9 in Fig. 
zon lAHl is seen to go to small 



1(b) 



1(b)) 



the internal hori- 
ler u (higher r). This is 



the manifestation of the antifocusing gravitational effects 
which are opposite to the effect 2), mentioned above in 
the discussion of the positive energy density pulse (see 



section rVII|) . 

The behavior of the mass function inside the pulse in 
the region 5 < < 9, (shown in Fig. 2(a) along lines 
of constant u), is also opposite to the case of an ingoing 
pulse with positive energy density described in section 
IVIIl in that the mass function decreases as the pulse 
enters the BH. This is also the case for pulses with higher 
(negative) energy contents (Fig. | 2(b)| . In Fig 



2(a) 



it 

can be seen that the mass function along the line u — Q 
(where we formulate the initial conditions) is lower than 
along the line w = 23 (slightly outside of the OAH). This 
is the manifestation of the scattering processes outside 
the BH. Because of the scattering process, only part of 
the energy along u = Q will actually reach the OAH near 
M = 23, the rest will be converted into outgoing radiation 
that will forever escape from the BH. 

The consequence of the concentration of the incoming 
energy near the inner Cauchy horizon, is still the same 
as it was in the case of the positive energy density pulse. 
This concentration can be seen in Fig. 3(a) where lines 
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FIG. 1: Positions of tlie apparent horizons for the case : A$ = 0; A<s, = 0.01, (v*o = 5, v*i = 9). When pulse of the negative 
energy crosses the outer apparent horizon (OAH), it becomes smaller and the inner apparent horizon (lAHl) becomes larger. 



1.005 
1.004 
1.003 
1.002 
1.001 
1 

0.999 
0.998 
0.997 
0.996 
0.995 



u=29- 
u=28- 



10 



u=27 




u=26 




u=25 


ction 


u=24 


Massfur 



u=0 



15 20 
V 



25 



30 



1.15 

1.125 
1.1 

1.075 
1.05 

1.025 
1 

0.975 
0.95 
0.925 




10 



15 



20 
V 



u=29 



-u=27 



— u=26 



u=25 



-u=23 
-u=0 



25 



30 



2(a): Case; = 0; A* = 0.01, (d*o = 5, 11*1 = 9). 



2(b): Case: Ag> = 0; = 0.045, (t;*o = 5, -y^i = 9). 



FIG. 2: The mass function along lines of constant u for two cases of weak pulses. One can sec the decrease of the mass function 
caused by the pulse of negative energy and the increase of it that results from the mass inflation. 
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3(a): Case: A<s, = 0; Aip = 0.01, (s;*o = 5, fwi = 9). 



3(b): Case: A$ = 0; A^s, = 0.045, (I'^o = 5, v^i = 9). 



FIG. 3: V versus r along lines of constant u for two cases of weak pulses. All lines with u < 24 comes to the same r when 
V — >■ oo. 
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4(a): Case: = 0; Aig = 0.01, (1)4,0 = 5, f*i = 9). 



4(b): Case: = 0; A^, = 0.045, (-Uspo = 5, -y*! = 9). 



FIG. 4: Kretschmann curvature scalar along lines of constant it for two cases of weak pulses. All lines inside the BH {u < 24) 
is seen to increase when v ^ 00. 



of r versus v along u — const asymptotically tend to 
r = rcauchy for u big enough. 

The sharp increase of the mass- function for the greater 
V inside the BH (for v > 12) in Fig. 2(a) is the result of 
the mass inflation which still works in the case of an 
exotic radiation pulse. 



The Kretschmann scalar along u = const is shown in 
Fig. 4(a) It is seen to increase with v for high values of 



u = const. The exponential increase of the Kretschmann 
scalar with u is a fingerprint of coming to the weak sin- 
gularity at u = 00. Thus in this case, probably the weak 
singularity exists at the Cauchy horizon. 
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5(a): Ingoing flux along lines of constant u. 5(b): Outgoing flux T^^ along lines of constant u. 

FIG. 5: In- and outgoing ^' fluxes of the negative energy for the case: A<s, = 0; Axjj = 0.01 (i;*o = 5, v^i — 9). The initial 
ingoing pulse can be clearly seen (fig. 5a, between 5 < w < 9), as can the scattered fluxes of the initial pulse. 



Fig. 5(a) represents the evolution of the in- flux, Tyy, of 



the negative scalar energy into the BH. Fig. 5(b) shows 
the corresponding T„„ out-flux which arises as a result of 
the flux being scattered by the space-time curvature. 
The tails of the Tw flux are the result of backscattering 
of the Tuu flux and so on and so forth. 
In both Fig. 



5 a 



and Fig. [^(b)] - 
nances arising as a result of the scattering process. These 
resonances are typical also for the cases of the irradiation 
by the more intense pulses. 

Now let us consider the case of a more intense pulse 
of negative energy, but still not strong enough to destroy 
the BH. We consider the case 

=0; A^ = 0.045, (i;*o = 5, = 9) (20) 



2(b) 



First of all we note the following; In Fig. 
see that at u = the value of the mass function after the 
pulse V > 9 corresponds to the external mass m « 0.945. 
This is less than q = 0.95 which means, that if this value 
does not change when one comes to the OAH, the BH 
must be destroyed. But because of the scattering process, 
part of the negative energy of the pulse will be scattered 
away during the propagation of the pulse from u = to 
the OAH and the amount of the negative energy pumped 
into the BH will be smaller. As a consequence, the flnal 
mass of the BH after the end of the whole process will be 
m w 0.959, and the position of the OAH process will be 
u w 23.6, see Fig. IHl Thus the BH survives. 

The reason for the increase of the mass function (in 
Fig. I 2(b)| in the region just after v ~ 10 for great u 
is mainly the additional flux of the negative energy r„„ 



which arise because of the scattering effect (analogous to 
the scattering demonstrated in Fig. [5]) and for bigger v 
it is a mild mass inflation. 



The properties of the Kretschmann scalar inside of the 
BH is now more complicated. At big u it shows growing 
oscillations with v, see Fig. 4(b)[ These oscillations are 
the consequence of the oscillations of the T^u and Tyy 
fluxes which are the analogues to those in Fig. [5j 



Finally in Fig. [5] one can see the general structure of 
the R- and T-regions. For all situations we will call the 
OAH the horizon which is practically the border of the 
BH. The lAHl is the border between the T_-region and 
the internal i?-region and its continuation. The horizon 
IAH2 is the border between the internal R and T+ regions 
and its continuation. All these horizons can be seen in 
Fig. [6l Most interesting is the fact that the border IAH2 
is not along a null-geodesic line (Cauchy Horizon) but 
comes to be visible inside the domain of integration. 



On the basis of our computations and their analysis we 
may guess at the Penrose diagram in Fig. [7] for the case 
when the BH survives after being irradiated by a pulse 
of the exotic scalar field. 
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FIG. 6: Lines of constant r (thin lines, decreasing from bot- 
tom to top), positions of R- and T-regions and positions of 
the apparent horizons (thick lines) for the case of Ag> = 0; 
A\s! = 0.045, (11*0 ~ 5, 11*1 =9). In this case, the BH sur- 
vives as can be seen by the OAH and lAHl not meeting. 




FIG. 7: The Penrose diagram for the case when the BH sur- 
vives after being irradiated by a pulse of the exotic scalar field. 
Here: 1 - strong, r = 0, timelike singularity, 2 - weak singu- 
larity, 3 - horizons, 4 - narrow pulse of the exotic scalar field, 
5 - light infinity, R' - our universe, R" - another universe, 6 
- the event horizon. 



B. Destruction of the BH 

Next we consider the following cases with (higher) neg- 
ative energies of the exotic pulse: 



A. A^ =0; A* = 0.05045 

B. A$ =0; A* = 0.0513 

C. A^ =0; A^= 0.0700 



(21) 
(22) 
(23) 



8(a)]r8(c)] show the 



In all cases v^o = 5, v^i = 9. Figs, 
further evolutions of the general picture of the R and T 
regions for these cases. 

In all these cases, the total power of the exotic pulse 
is large enough to reduce the mass of the object in 
the critical region it w [23.2; 23.9] ( where the OAH 
formed for the previous cases) , to below the critical value 
merit = q = 0.95, see Fig. [SI Thus even after the reduc- 
tion of the power of the pulse during the propagation 
from u — to u 23.2 due to the scattering, it is strong 
enough to destroy the BH. 

This means that the outer and inner apparent hori- 
zons should meet and disappear. This process is seen 
in Fig. 



a) and even more clearly in Fig. | 8(b)[ which 



shows the slightly stronger case 

For these cases, for high v, the BH is converted into an 
object where the outer R and inner i?-regions are con- 
nected. Now the test photons u — const for all value of 
u go to bigger r, when w — > cxd. 

For all cases (j2m23|) the mass-function does not 
demonstrate a mass inflation effect, see Fig. [HI We note 
that for the cases (I^Tj) and the mass function for big 
V and big u does becomes greater than m = 0.95. How- 
ever, this is related to the scattering of the exotic scalar 
field outside to bigger r. Of cause this is possible in the 
case of dynamical BHs. In this region (big v and big u) 
we are at big r, definitely outside the (dynamical) BH 
which is at smaller v. 

For the cases (|2m22p the Kretschmann scalar does not 
increase with u, but rather oscillates (see Figs. fTU]) . These 
oscillations are the consequences of the oscillations of the 
r„.u-in flux (and also Tti„-out flux), which we described 
for the pulses with the smaller power, see Fig. 5(a)|| 5(b)" 



The fact that the Kretschmann scalar does not increase 
exponentially with v, together with the fact that the mass 
function does not increase exponentially, indicates the 
absence of space-time singularities (excluding the time- 
like r = singularity of the Reissner-Nordstrom BH be- 
yond the top- left corner of the computational domain) . 
Finally, in Fig. 8(c) is shown the horizon structure for 



the case (^5)). 

In this case, the exotic pulse is so powerful that the 
inner and outer apparent horizons meet quite fast and 
after this, there are not any horizons or borders left at 
all (to the right of u ~ 7.5). The outgoing signal u — 
const can freely propagate away through the i?-region. 
Also, the mass function becomes essentially smaller then 
q (see Fig. 9(c) |, while the Kretschmann scalar becomes 



smaller and smaller with bigger v (see Fig. [TUT c)). There 
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8(a); Case: A* = 0; Aq, = 0.05045. 



9(a): A* = 0.05045 
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8(b): Case: A„ = 0; A,, = 0.0513. 



9(b): A* = 0.0513 
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8(c): Case: A* = 0; A,p = 0.070. 



9(c): A* = 0.070 



FIG. 8: Lines of constant r (thin lines, decreasing from bot- 
tom to top) , positions of R- and T- regions and positions of the 
apparent horizons (thick hncs) for the cases a) Ais = 0.05045; 
b) A* = 0.0513; c) A* = 0.070. In aU cases (t;*o = 5, 
Vijfi =9) and in all cases the BH is destroyed as can be clearly 
seen by the meeting of the horizons. 



FIG. 9: The mass function along lines of constant u for the 

cases: a) A* = 0.05045; b) A* = 0.0513; c) A* = 0.070. In 
all cases A$ = 0; (w*o = 5, v*i = 9) and in all cases the BH 
is destroyed. 
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10(a); Case: A^g = 0.05045 




10(b): Case: Aq, = 0.0513 




10(c): Case: A^ = 0.070 

FIG. 10: Kretschmann scalar along lines of constant u for the 
cases: a) A* = 0.05045; b) A* = 0.0513; c) Aq, = 0.070. In 
all cases A^ — 0, {v^o = 5, = 9) and in all cases the BH 
is destroyed. 




FIG. 11: The Penrose diagram for the case when the BH is 
destroyed after being irradiated by a pulse of the exotic scalar 
field. Here: 1 - strong singularity r = 0, 2 - horizons, 3 - 
narrow pulse of the exotic scalar field, 4 - light infinity, R' - 
our universe, R" - another universe. 



are not any singularities present ( except for the timelike 
r = singularity). 

Fig. [TT] represents the Penrose diagram that is con- 
firmed by our numerical simulations for the cases when 
the BH is destroyed by the radiation. 

At the end of this section we note the following. When 
the Reissner-Nordstrom BH is irradiated by a pulse of 
the exotic scalar radiation, the OAH becomes smaller 
(or disappears completely) and part of the outgoing ra- 
diation from the T_ region can go to the outer i?'-region 
in our Universe. This radiation may come into the T_ 
region from the i?"-region that belongs to another Uni- 
verse, which is the counterpart of the outer i?'-region of 
Fig.[Tl]in our Universe (from the left hand side of Fig. [TT] 
outside the computational domain). This means that it 
is possible for some radiation from the other Universe to 
come to our i?'-region. The propagation of the radiation 
in the opposite direction, from our i?'-region to the i?"- 
region in the other Universe, is still impossible. We call 
such an object a semitraversable wormhole (soj . 



IX. THE CASE OF THE IRRADIATION BY 
BOTH NORMAL AND EXOTIC PULSES 



Finally we consider the case: 



0.4, 

:0.2, 



5, 
10, 



7) 



12). 



(24) 
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FIG. 12: Lines of constant r (thin lines, decreasing from bot- 
tom to top), positions oi r = 0, K = Kpianck, apparent hori- 
zons and R- and T-regions for case: j4$ — 0.4, {v<s>o ~ 5, 
= 7); Aij! — 0.2, (u^o = 10, u*! = 12) in which the BH is 
destroyed. 



In this case, the Rcissner-Nordstrom BH is irradiated 
at the beginning by a strong pulse of the normal scalar 
field which causes a strong spacelike r — singularity to 
arise in the T_ -region inside the BH (similar to the cases 
investigated in [33l). After that, at higher v, the BH is 
then irradiated by a strong pulse of the exotic scalar field. 

For this case, one can see strong nonlinear processes. 
Fig. [12] shows the general structure of the apparent 
horizons and the position of the (spacelike) singularity 
r — 0. There is also a line showing the position of 
the Kretschmann scalar equal to its Planckian value of 
K = {Gh/c^)-'^ w 10^3ic„i-4 Beyond this border the 
classical General Relativity is not applicable and we con- 
sider this region as a singularity from the classical point 
of view. Our calculations are thus not reliable for this 
region and we will not consider it (upper-right part of 
Fig. [121 above the line K = Kpianck)- We would like to 
do the following remark about the iiT-scalar. This scalar 
has dimension cm~'^ in the CGS -units system. In the 
paper [33] and in this paper we represent all figures in di- 
mensionless units. For the figures, we draw the i^T-scalar 
(and the value of the Planckian ivT-scalar) for the case 
when the linear scale (the size of the BH in this paper) 
is equal to unity. 

The Kretschmann scalar is represented in Fig. [T31 For 
large u > 14, the Kretschmann scalar grows catastroph- 
ically fast and reach the singularity K — Kpianck- This 
is the consequence of the influence of the first pulse, with 
the positive energy density. This is the same process 
which we analyzed in [s^, see e.g. Fig. 13 in [ssj. Along 



FIG. 13: Kretschmann scalar versus v along lines of constant 
u for case: A<s, = 0.4, {v-s,o = 5, «*i — 7);^* = 0.2 (i;*o ~ 10, 
= 12). Also, the location of the Kretschmann singularity 
Kpianck is plotted. 



u — 12, the Kretchmann scalar decreases with increas- 
ing V, this is of course as expected as this line is always 
outside of any horizon (fig. [T^ . The behavior along the 
line u = 13 is somewhat in between, for small v {v ^ 14), 
the Kretschmann scalar shows a slight tendency to in- 
crease. However, for larger v, the BH is destroyed and the 
Kretchmann scalar shows a decreasing behaviour. Also 
cf. with the discussion of the mass function below. 

In Fig. [Til (showing r versus v), the test photons along 
u — 18 reaches r = 0. Photons in the range 13 < u < 17 
come close to r = and after that, they propagate to 
larger r (after the meeting with the pulse of the exotic 
scalar field atl0<t;<12). However, it is reminded that 
photons along u > 14 meet with K = Kpianck and thus 
effectively reach the singularity before they can escape 
to higher r. The photons along m = 12 constantly go to 
bigger r because they are always outside of any horizon. 

In Fig. [ini it is seen that the mass function goes up 
for all lines of u = constant between 5 < f < 7, this 
is merely the trivial effect of mass being pumped into 
the BH by the positive energy pulse. In addition to this, 
the mass function for m > 14 continues to increase even 
after the passage of the positive energy pulse, due to the 
processes described in section IVIII until they reach the 
K = Kpianck singularity. 

As the negative energy pulse goes into the BH, the 
mass function along lines u < 13 decreases due to the 
negative energy contents of the second pulse. This sec- 
ond pulse, causes the mass function to become less than 
q in the region near the outer horizon and it leads to the 
meeting of the inner and outer apparent horizons and the 
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FIG. 14: r versus v along constant lines of constant u for case: 
Ai, = 0.4, (uso = 5, usi = 7); Aq, = 0.2, (v-ao = 10, u*i = 
12). All lines for u > 17 can in principle escape to infinity, 
except that most of them hit the Kretschmann singularity 
(see text for details). 

disappearance of the BH around v ~ li (see figure [T2l) . 
Along u = 12, the mass function continues to decrease be- 
cause of the scattering processes of the exotic scalar field. 
Along u — 13, however, we see that the mass function in- 
creases drastically near w w 12. This growth is probably 
caused by the mass inflation effect, since in this region 
we are close to the border of the BH, but still inside it. 
The apparent horizons do not meet and the BH is not de- 
stroyed until around v ~ 14. In the same region, one can 
also see a modest growth of the Kretschmann scalar in 
Fig. [T31 However, as the BH is destroyed around v « 14, 
we see that the mass function halts its growth and the 
Kretschmann scalar starts to decrease. 



X. CONCLUSIONS 

The processes arising when a Reissner-Nordstrom BH 
is irradiated by a pulse of an exotic scalar field with a neg- 
ative energy density have been analyzed. We have per- 
formed the corresponding numerical computations using 
a numerical code specially designed for the purpose. It 
was demonstrated that these processes are quite different 
from the processes arising in the case of the irradiation of 
a Reissner-Nordstrom BH by a pulse of a normal scalar 
field. 

The evolution of the mass function and the 
Kretschmann scalar demonstrate that in the case of the 
exotic scalar field, the evolution does not lead to the ori- 
gin of a strong spacelike singularity r = in the T-region 



FIG. 15: Mass function along lines of constant u for case: 
Ai, = 0.4, (u<5o = 5, = 7); Aq, = 0.2, («so = 10, vq/i = 
12). Effects of both the normal and the exotic scalar pulse 
are clearly seen. 



as was the seen in the case of irradiation by the normal 
scalar field I3il] . 

The numerical calculations demonstrate the manifes- 
tation of the antifocusing effects in the gravity field of an 
exotic scalar field with a negative energy density. 

When the power of the exotic pulse with negative en- 
ergy density is great enough, the mass function becomes 
less than the charge q near the outer apparent horizon. 
As a result the BH disappears. This process was analyzed 
in detail. 
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Appendix A: Convergence of the code 

In this appendix, we demonstrate that our numerical 
code is converging (to a physical solution) when includ- 
ing the effects of nontrivial $ and ^ fields. The results 
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presented here are very similar to those presented in [3J| 
as could be expected, since the codes in that paper and 
this are essentially identical. For this reason, for further 
details of the code, we refer to [33 |. 

The initial conditions for the tests in this Appendix are 
similar to those described in section IIXI i.e. the initial 
conditions are those of a Reissner-Nordstrom BH, per- 
turbed by first an in- falling $ field (from v — [5; 7] with 
amplitude: A$ = 0.4) and subsequently by an in-falling 

field (from i; = [10; 12] with amplitude: = 0.2). 
Our computational domain for the convergence tests are, 
as for all simulations in this paper, in the range v = [5; 30] 
and u = [0;30]. This setup is the most complicated one 
that we have done in this paper (in that it incorporates 
the dynamic effects of both <i> and ^ fields) and leaves 
no trivial terms left in the evolution equations, thus it is 
a good test for the convergence of our code. 

Here we demonstrate the convergence of the code by 
comparing a series of simulations with varying numerical 
resolution. For the tests in this Appendix, we do a to- 
tal of 6 simulations, with each simulation changing the 
AMR base resolution and refinement criteria in such a 
way as to mimic doubling the numerical resolution from 
one simulation to the next (for further details and more 
thorough discussions, see |34l i44|). 

To limit the number of plots in this Appendix, we con- 
centrate on displaying convergence results along the line 
u = n (however, it should be noted that a large number 
of convergence tests has been carried out and that the 
results presented in this Appendix are representative of 
the convergence properties of the code in other parts of 
the computational domain and for other configurations 
of the initial data) . As can be seen in figure [TH this line 
is very close to the r = singularity (where we can ex- 
pect strong dynamics) as well as inside the Kretschmann 
singularity and thus convergence along this line is a good 
indication of the convergence behavior throughout the 
computational domain. 

Along this line, we calculate the relative convergence 
between two simulations (one with a numerical resolution 
twice that of the other) relative to a simulation with very 
high resolution: 



''2N\ 



'^HighRes I 



(Al) 



where x]^ denotes the dynamic variable x at the i-th grid 
point of simulation with resolution N and where a^jj jg/i7?es 
denotes the dynamic variable of the same i'th point for 
a simulation with the highest numerical resolution done 
by us. Obviously, this expression only makes sense for 
those i points that coincide in all simulations. 

Figures 16(a)|fl-6(d) show the relative convergence, ^, 
for the dynamic variables r, a, <i> and 5* respectively. The 
lines in the figures are marked by their numerical reso- 
lution measured in terms of the most coarse resolution 
Nq, the high resolution simulation used to calculate ex- 
pression (|Aip , has a numerical resolution of 32 times the 



base resolution, i.e. N = 32No. 

From these figures, it is clearly seen that the four dy- 
namic variables are converging for simulations of increas- 
ing resolution. Furthermore, since we plot the relative 
convergence of the dynamic variables, we see that the 
relative change between the two highest resolution simu- 
lations show that the variables change 0.1% or less, which 
must be considered a quite acceptable convergence. We 



note that a closer analysis of the data in figures 16(a 



16(d) has revealed that the dynamic variables are indeed 



converging with second order accuracy as was expected 
based on analyses and tests performed in [3^ [s^j . 

However, of course it is not enough to demonstrate 
that the simulations are converging, they must also con- 
verge to a physical solution, i.e. the residuals of the con- 
straint equations (eqs. ^ and ^) must converge to 
zero. To demonstrate this, we calculate the relative con- 
vergence of the constraint equation residuals, (relative to 
the Einstein-tensor, eqs. (|2a|) and ([2b)) respectively) in a 
similar way to eq. (jAll) : 



N) 



N\ 



\G' 



(A2) 



HighRes I 



where denotes the residual of the constraint equation, 
{Cuu or C„„), at the i'th point for simulation with reso- 
lution N and where G^^jg^^g^ denotes the corresponding 
Einstein-tensor component {Guu or Gw respectively) at 
the same point. 

The relative convergence of the residuals of the con- 
straint equations are demonstrated in figures 16(e) and 
16(f) where it is seen that they converge towards zero 



for higher resolution simulations. This indicates that not 
only are the numerical solutions converging for simula- 
tions of higher resolution, but that they are indeed con- 
verging towards a physical solution. The convergence 
plots for the constraint equations exhibit significant more 
numerical noise than the convergence plots for the basic 
dynamic variables. However, it is reminded that the con- 
straint equations are derived quantities of the basic dy- 
namic variables, therefore the constraint equations will 
be less accurate and include more numerical errors such 
as truncation errors from calculating the derivatives of 
the basic dynamic variables and errors associated with 
calculating derivatives in the AMR mesh. Nevertheless, 
figs. 



16(e) and 16(f) 



indicate that the solutions are 
indeed converging to the physical solutions for high res- 
olutions. 

Finally it should be noted that the convergence results 
presented in this appendix are not the only convergence 
tests that we have performed, they merely represent typ- 
ical results of the convergence behavior of the code. For 
all results presented in this paper, we have performed a 
large number of simulations with varying resolutions to 
ensure that the results had converged to their physical 
solution. 
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constraint equations along line of u = 17, see text for details. 
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